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Molecular  Oroital  Study  of  Acyclical  Nitramine  Compounos 
I.  INTRODUCTION 

Nitramines  (nitro-amines)  are  chemical  compounds  which  have  the  generic 
structure  shown  in  Fig.  1.  There  is  a  nitro  group,  NO2,  double  bonded  to 
a  nitrogen  atom,  wnile  the  moieties  X  and  Y,  which  may  be  members  of  a  ring 
or  radicals,  complete  the  amine  structure;  a  resonance  structure  for  the  N-0 
bonds  is  understood. 

Such  compounos  are  of  both  theoretical  and  experimental  interest  because 
a  number  of  them  are  explosive,  while  others  may  serve  as  models  of 
explosives.  Figure  2  contains  the  structures  of  the  explosive  RDX,  and  the 
model  compound  uimethylnitramine  (OMNA).  The  1,  2,  6  positions  of  RDX, 
shown  above  the  dashed  line,  are  modeleo  by  the  part  of  DMNA  above  that 
line.  The  moieties  X  and  Y  correspond  to  the  CH2  carbenes  in  both  RDX  and 
DMNA.  In  the  former,  they  are  members  of  the  ring;  in  the  latter,  they 
combine  with  the  remaining  H  atoms  to  form  the  CH^  methyl  radicals. 

'^/Because  the  reactions  of  such  molecules  are  complicateo,  it  may  happen 
that  the  first  steps  of  their  oecompositions  are  not  amenable  to 
experimental  study,  and,  in  order  to  attempt  to  obtain  an  understanding  of 
the  complete  dynamics,  it  is  necessary  to  perform  theoretical  calculations. 
The  small  value  of  the  ratio  of  electron  to  proton  masses,  1/1823,  leads  to 
a  decoupling  of  the  electronic  motion  from  that  of  the  heavy  particles,  the 
Born-Uppenheimer  approximation.  The  quantum  mechanical  treatement  of  the 
electronic  motion  for  fixeo  positions  of  the  nuclei  forms  the  basis  of 
quantum  chemistry  which  generates  a  potential  energy  surface  that  determines 
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tne  heavy  particle  motion  classically,  semi-classical ly,  or  quantum 
mechanical ly» 


in  Sec  II,  the  Hartree-Fock  (HF)  self-consistent  field  (SCF)  method  is 
briefly  oescribed.  It  underlies,  in  some  fashion,  all  forms  of  quantum 
chemistry.  The  reduction  of  the  HF  equations  by  means  of- atomic  orbitals 
(AO's)  is  described  as  is  the  proolem  of  multifolo  proliferation  of  matrix 
elements.  In  Sec, III,  the  reouction  of  the  full  SCF  Hamiltonian  to  the 
.'iliiOO/S  ana  forms  is  oiscussed,  ana  the  empirical 

parameterization  of  matrix  elements  is  described.  In  Sec.  IV,  the  results 
obtained  for  the  ground  state  properties  of  several  acyclical  nitramine 
compounds  using  the  MINDO/3  and  MNDO  Hamiltonians  are  compared  with  each 
other  and  with  experiment.  Conclusions  are  contained  in  Sec.  V. 


II.  THE  SELF-CONSISTENT  FIELD-LINEAR  COMBINATION  OF  ATOMIC 
ORBITALi-MOLECULAR  ORBITAL  (SCF-LCAO-MO)  METHOD 


The  complete  non-relativistic  molecular  Hamiltonian  may  oe  written  as 


H=i;  T+e  „  V  +e  t  +e 
a  a  a>S  aS  1  i  I, a 


V  +  E 

i  a 


V 

ij  ’ 


(1) 


where  the  sum  on  Greek  indices  goes  over  the  atomic  nuclei  and  that  on  Latin 
indices  goes  over  the  elctrons.  The  quantities  T  are  the  kinetic  energies 
and  the  symbols  v  denote  the  Coulomb  potential  energies.  In  terms  of 
natural  units  such  that 


£q  (1)  acquires  the  form 


H=(l/MjZ(p2/2A  )  +  Z  Z-|r  -R„r^ 
A  a  a  a  a>8  a  8  a  6 


+Z(p  2/2)-. E  Z  IttR  +  ,E.  |r.-r. 

i,aa'  xa'  i>j  '  i  j '  ♦ 


±  I  “1 


(1.1) 


where  the  quantities  p(P)  are  the  momenta,  the  quantities  z  the- atomic 
numbers,  tne  nuclear  masses  in  atomic  mass  units  are  A  ,  while  R  ana  r 
denote  the  positions  of  the  corresponding  particles.  The  occurrence  of  the 
factor  =  (1/1823),  the  inverse  of  the  atomic  mass  unit  expressed  in 

natural  units,  suggests  the  Born-Oppenheimer  approximation,  the 
consideration  of  Eq.  (1.1)  for  fixed  positions  of  the  nuclei.  Treating  the 
nuclear  coordinates  as  parameters  leads  to  the  fundamental  equation  of 
quantum  chemistry 

H^(?,R)f  (?,?)=u  (5)¥  (?,R), 
e  m  mm 

Where  rand  1  denote  all  of  the  and  ,  respectively  ana 
H^(;,R)-£(V^^/2+EZ^|?j-S„r‘) 

Js  V6lV®6!"‘  .  (3) 

The  functions  which  diagonalize  Eq.  (2)  are  adiabatic  basis 

functions,  and  the  eigenvalues  are  adiabatic  electronic  energies.  In 

order  to  solve  the  complete  problem  implied  by  Eq.  (1),  non-adiabatic 
correction  terms  must  be  added  to  the  potential  u  in  order  to  obtain  a 
solution  for  the  nuclear  part  of  the  total  wave  function.  Alternatively, 
one  may  choose  to  define  electronic  energies  as  the  eigenvalues  of  a 
Hamiltonian  other  than  that  given  by  Eq.  (3)  in  which  case  the  functions 
corresponding  to  the  are  known  as  diabatic  basis  functions,  ana  the 
eigenvalues  as  diabatic  potential  energy  surfaces. 


m 


The  variational  method  supplies  a  means  of  constructing  approximate 
solutions  of  Eq.  (2)  according  to  the  prescription  that  the  energy  functional 


be  stationary  with  respect  to  variations  of  the  functions  |'"^>(or<fml) 
suDject  to  the  subsidiariary  condition 


1  . 


The  variational  method  becomes  the  HF  method  when  the  are  restricted 

to  be  single  Slater  determinants  composea  of  one  electron  spin-orbital 
functions.  For  the  ground  state  of  a  closed  shell  system  of  N  electrons,  it 
is  convenient  to  suppress  the  subscript  m  and  write  |4'ni>  in  the  form 

|f>  =  (N!)"^ 


where  the  ni  assume  the  values  or  according  to  the  z-component  of  the 
spin  of  the  electron.  The  |(()i>  are  a  set  of  single  particle  orbitals, 
and  the  sum  extends  over  the  permutations  P  of  parity  (-)^  that  may  be 
formed  from  the  product  given  in  a  standard  order.  Because  the  resulting 
equations  are  non-linear,  they  are  solved  by  interation  until  convergence  of 
two  (or  more)  consecutive  iterations  to  within  a  given  tolerance  is 
achieved,  i.e. ,  the  solution  is  self-consistent  and  the  resulting  energy  is 
said  to  have  been  calculated  by  the  HF-bCF  method.  The  method  becomes  the 
HF-SCF-LCAU-MO  method  if  a  linear  combination  of  atomic  orbitals  (LCAOj 
centered  on  the  nuclei  a  are  chosen  for  the  molecular  orbitals  (hO)  j((>j^>in 


which  case 


The  inaex  p=(:i,a)  where  X  labels  the  atomic  orbitals  centerea  on  the 
nucleus  ct ,  ana  the  \xp>  are  the  normalized  atomic  orbitals.  The  upper 
limit  to  the  sum,  K,  has  a  value  not  less  tnan  N/2.  if  it  were 
allowea  to  be  infinite,  the  |c|)i>  woula  be  a  complete  set,  ana  the  solution 
of  Eq.  (4)  would  be  at  the  hF  limit.  If  k  were  to  correspond  to  the  minimal 
number  of  functions  neeaed  to  describe  the  valence  electrons.  Is  for  H  ana 
2s,  2p^,  2py  ana  2pz  for  first  row  atoms,  the  basis  set  is  minimal.  A 
minimal  basis  set  is  employed  in  the  calculations  to  be  discussed. 


The  definition  (6)  does  not  imply  that  the  |>}>^>  form  an  orthonormal  set, 
but  a  unitary  transformation  enables  them  to  be  written  as  such.  The 
variation  of  the  anti-symmetrized  product  may  be  written  as  a  variation  of 
the  Cp^.  The  introduction  of  Lagrange  multipliers  in  order  to  enforce  the 
requirement  of  orthonormality  of  the  |<()^>,  ana  variation  of  the  leaa  to 
the  HF-SCF-LCAO-MU  equations 

i  C,  -  cM  ,  l^i^N/2  (7) 

i  i  i 

where  c  is  a  column  vector  with  K  components 


and  S  is  the  overlap  matrix  with  elements 
s  =  <x  I  X  > 

pq  ^p  ' 

The  Fock  matrix,  ?,  has  elements 


"pq  “  ' 


(9) 

(10) 


where 


V,/2-  Z  Z  /|r,-R 
1  a  a  '  1  a 


+  2  Z„/|r  -R„1 

a>3  a  S  a  6 


(") 


the  Coulomb  operator  is 


-1 


j.  =  <6^(2)|r^^iA.(2)>  , 


(12) 


ana  the  exchange  operator  is  definea  by 


t^ilXq>  =  ir^2lXqC2)>Ui>  • 


(13) 


The  quantities  ei  are  the  molecular  orbital  single  particle  energies.  The 
factor  2  arises  in  tq.  (lOJ  because  all  states  are  doubly  occupied  in  the 
closed  shell  system.  Because  the  matrix  elements  of  Eqs.  (12)  and  (13) 
involve  matrix  elements  over  molecular  orbitals.  Equations  (7)  are,  in  fact, 
non-linear  and  of  the  third  degree  in  the  expansion  coefficients  c  .. 


Equations  U)  have  a  deceptively  simple  appearance.  Not  only  has  the 
non-linearity  of  the  problem  been  surpressed  but  also  there  is  a  vast  number 
of  two  electron  integrals  to  be  evaluated.  Neglecting  possible  reductions 
in  the  number  of  integrals  due  to  molecular  synmetry,  there  are 
Nj,  =  (K/4)(K+1)  [(K/2)(K+l)+l] 

distinct  two-electron  integrals  for  a  K-dimensional  basis  set  calculation  of 


the  F  .  In  the  case  of  a  molecule  such  as  DMNA,  a  minimal  basis  set  for 
pq 


only  the  valence  electrons  has  30  members  and  N2q=1083451  Many  of  these 
may  be  expected  to  be  small,  however,  and  this  leads  to  the  introduction  of 
semi-empirical  approximations. 


III.  SEMI -EMPIRICAL  HF-SCF-LCAO-MO  APPROXIMATIONS 


Reducing  the  number  of  two-electron  integrals  that  require  evaluation 
simplifies  the  computation  of  the  solutions  of  Eqs.  (7).  The  most 
optimistic  assumption  about  the  matrix  element 


=  (pq|rs) 


(14) 


6 


IS  that  if  either  or  both  of  (p,q)  and  (r,s)  refer  to  different  atomic 
orbitals,  oscillations  in  the  integrand  due  to  differing  nodal  structures 
will  cause  the  matrix  element  to  be  small.  The  zero  differential  overlap 
approximation  then  sets 


(pq|rs)  =  (pp|rr)6  6 

pq  ri» 


=  y  .6  6  .  (p,q)  Ca;  (r,s)  C2 

ctS  pq  rs 

Further  approximations  must  be  introduced  in  order  to  construct  a  theory 
invariant  to  unitary  transformations  of  both  the  MO's  and  AC's 


Two  particular  sets  of  approximations  lead  to  the  construction  of  the 

2  3 

MINDO/3  and  MNDO  Hamiltonians,  which  however,  have  several  common 

features.  Both  of  them  incorporate  the  assumption  that  a  frozen  core  of 

inner  shell  electrons  completely  shields  the  bare  nuclear  charge. 

Similarly,  only  a  minimal  basis  set  of  AO's  is  retained  for  the  valence 

electrons;  thus,  hydrogen  is  represented  by  a  single  is  orbital,  and  the 

remainder  of  the  first  row  atoms  are  represented  by  the  four  orbitals  2s, 

2p  ,  2p  and  2p  .  The  orbital  exponents  are  adjusted  from  their  usual 
X  y  z 

values  by  treating  them  as  parameters  during  the  empirical  determination  of 
(pp|rr).  Both  models  neglect  all  three  and  four  center  integrals,  and 
presume  that  overlap  matrix  elements  enter  the  wave  function  normalization 
matrix  S  as  the  unit  matrix,  but  may  differ  from  the  unit  matrix  elsewhere. 
Monatomic  differential  overlap  is  retained  for  all  one  center  integrals. 


In  i'ilNOO/S*^  one-center,  one-electron  matrix  elements  are  related  to 
observed  ionization  potentials  and  spectroscopic  data,  as  are  the 
one-center  two-electron  repulsion  integrals.  The  attractive  integrals 


between  an  electron  centerea  on  atom  a  ana  the  core  of  a  different  atom  3 
are  assumed  to  depend  only  upon  the  pair  of  atoms  involved,  and  are 
uetermineo  by  a  phenomenological  form  wnich  reduces  to  the  classical  Coulomb 
atrraction  when  R,  ^  ,  the  distance  between  atoms  o.  ana  3,  becomes  large. 

Tne  two-center  core  resonance  integrals  are  assumed  to  be  given  by 


s 

pr 


(1^+13) f(R 


U6) 


where  S  is  the  overlap  between  orbitals  p  ana  r  centerea  on  atoms  ct 
pr 

ana  3,  respectively,  I  is  the  ionization  potential  of  atom  a(6),  and  f 

ctvp; 

is  a  function  to  be  determined;  in  practice  f  is  chosen  to  be  a  constant. 


B^g.  The  two-center  two-electron  repulsion  integrals  are  required  to  be 

determined  only  by  the  atom  pair  a  and  6  as  a  result  of  rotational 

5  -1 

invariance  ,  and  are  given  by  a  form  which  tends  to  R^g  for  large 

R^g,  and  on  the  average  of  one-center  two-electron  integrals  for  atoms  a 

ana  3  when  Finally,  the  core-core  repulsions  are  determined  by 

empirically  fitting  the  parameters  of  an  analytic  formula. 


3 

The  MNuo  Hamiltonian  extends  ana  generalizes  the  MlNDU/3 
Hamiltonian.  >iatrix  elements  of  the  form 

(pqirs)  (p,q)Ca,  (r,s)0  ^ 

are  retained  for  p#q,  r#s,  and  matrix  elements  of  the  form  (pp|rr)  are 

allowed  to  depend  not  only  on  a  ana  3  but  also  upon  the  values  of  p  and  r. 

(Nevertheless,  matrix  elements  of  the  form  (  pr|qs)  are  neglected  as  before. 

Tne  one-center  matrix  elements  are  treated  as  in  MINDU/3.  The  attractive 

integrals  between  an  electron  centered  on  atom  a  and  the  core  of  another 

atom  3  are  calculated  by  assigning  the  quantum  mecnanical  monopole 


TABLE  6.  Comparison  of  MNDO  results  with  experimental  data 
for  C^'(CH3  j  iN(N02) . 


Observable 

Experiment^ 

MNDO 

RvNN)(rt) 

1  .A69±O.CC5 

1.434 

R(iNO)  (Aji^ 

1 .209±0.002 

1.206 

R(CN)(AJ 

1.478±0.005 

1.492 

K^Ci.NJ(A) 

1.720±0.004 

1.735 

£.CNN(') 

112.9±1.5 

117.8 

0 

bCZNi'H  ) 

108.4±1.3 

111.0 

ICSLNC(°) 

115.0±1.2 

110.5 

bGhO(  ) 

128.5±0.9 

123.7 

lCZNH0(°) 

24.4±1.0 

50.7 

0 

bChNU(  ) 

27.C±1.0 

5.7 

AH(kcal  m“l) 

41.4 

X2 

128 

aRef.  9. 


^Averaqe  value. 


TABLE  5.  Comparison  of  MNDO  and  P'iINOO/3  results  with  experimental  data 
for  (CH2C^)(CH3)N(N02). 


Observable 

Experiment^ 

MNDO 

MINDO/3 

R(rt)b 

1 .426±0.010 

1.452 

1.431 

R(NO)(X)C 

1  .220±0.003 

1.208 

1.238 

R(CCil)lA) 

1.790±0.007 

1.814 

1.940 

Z.OiNO(  ) 

123.5±0.5 

123.3 

123.4 

Z.C^’NiM(°) 

116.2±0.5 

118.2 

117.5 

acnc(“) 

127.5±0.5 

117.8 

109.1 

tNCCi(°) 

112.3±0.8 

112.2 

120.5 

Z.C5nnU(°) 

35.0±1.0 

-0.3 

0.3 

0±1.0 

22.1 

-0.01 

AH(kcal  m"1 ) 

13.6 

127.8 

240 

398 

SRef.  9. 

t3R=(l/3)  [R(c5n)+R(C9n)+R(NN)]. 
^Averaije  value. 

^Discussed  in  text. 
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TABLE  4.  Comparison 

of  MNDO  ana  MINOO/3 

with  experimental  data  for 

(CH3)N(I'(02)2  .  a  non- 

-planar  molecule. 

Observable 

Experiment^ 

MNDO 

MINDO/3 

R(iNN)(A)a 

1  .480+0.005 

1.437 

1.355 

k(i\0)  (A)^ 

1.23U0.003 

1.203 

1.219 

R(CN)(A) 

1.494±0.006 

1.502 

1.460 

LMNN(^) 

117.0±1.1 

113.2 

122.7 

z.ONa('^)a 

132.0±1.0 

125.4 

132.3 

Z.CNN(°)3 

107.5±1.0 

114.6 

121.4 

^04n1n2,'^3(°) 

0.1±2.0 

-29.8 

67.8 

Z,o5m1o4.\)2(°) 

±180.0±1.0 

174.1 

-178.6 

Z.06n3n2i^1(°) 

0.0±2.0 

113.2 

19.0 

Z.o7n3o1n2(^) 

±180.0±2.0 

-179.0 

-177.3 

Z.CN2iN3n1(°) 

42.2±1 .0 

36.2 

5.3 

AH(kcal  m“l) 

53.9 

-36.4 

X2 

342 

318 

^Average  value. 

^All  experimental  data  are  from  kef.  15 

• 
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TABLE  3.  Comparison 

of  MNDO  ana  MINljO/3 

results 

with  experimental  data 

for  H(CH3)i^(N02)  . 

Observable 

Experimenta 

MNDU 

MI NOG/ 3 

R(NN)(A) 

1.381±0.006 

1.403 

1.297 

R(N0)tA)b 

1.228*0.003 

1.209 

1.230 

R(CN)(A) 

1.452*0.006 

1.479 

1.423 

Z.CNN(°) 

109.0-1.3 

119.6 

131.8 

^0N0(°) 

125.3*1.0 

123.7 

131.5 

^CNH(°) 

119.0*1.0 

110.3 

109.5 

/.CNN0(°) 

28*1 

35.8 

0.052 

AH(kcal  m"1) 

18.8 

-25.1 

X2 

39.9 

205.0 

aRef.  14. 

^Average  value. 
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TABLE  2.  Comparison  of  MNOO  and  i'^INDG/3  results  with  experimental  data  for 
(CH3)2N1N)2). 


Observaole 

Experiment^ 

MiMOO 

MINDC/3 

R(NN)(S) 

1.382*0.003 

1.395 

1.311 

R(N0)(]i;b 

1.223*0.002 

1.209 

1.232 

R(CR)(A)b 

1  .460*0.003 

1.486 

1.446 

R(CHj(^)b 

1.121*0.005 

1.115 

1.118 

z.QN0(°) 

130.4*1.5 

122.8 

130.2 

116.2^0.3 

117.9 

122.0 

aC5n2n1o3(°) 

0.0±1.0b 

19.0 

1.10 

Z,C6n2n1o3(°) 

180. 0*1. oa 

178.9 

180.0 

Z,u4n1n2o3(°) 

180.0*1 .Ob 

175.8 

180.0 

ii(O) 

4.61*0.066 

4.51 

4.19 

A H(kcal  m"l ) 

-16.9*1.0^ 

22.9 

-9.6 

X2 

210 

96 

3A11  experimental 

data  are  from  Ref. 

9  unless  otherwise 

notea. 

•^Averayed  over  all  bonds. 

CRef.  10. 

^Experimental  uncertainty  assigned  in  the  present  work. 
eRef.  11. 

^Ref.  12. 


TABLE  1.  Comparison  of  MNOO  and  MINDO/3  results  with  experimental 
data  for  H2N{N02). 


Observable 

Experiment^ 

MNDO 

MIhDO/3 

R(NhKA) 

1.427±0.002b 

1.406 

1.287 

R(N0)(A) 

1. 206^0. Oic 

1 .209*3 

1.230 

R(NH)(A) 

1.005±0.010 

1.018*3 

1.027*3 

115.2±2.0 

107.4 

108.3 

130.1±0.25 

124.1 

132.2 

0 

) 

51.8±1.0 

45.1 

0.03 

0  . 

Ln  ) 

0±0a 

5.8 

0.00 

/iH(kcal  m-1) 

-13.7±1.0e 

19.4 

-29.8 

u(D) 

3.57±0.05 

3.84 

3.95 

X2 

238 

998 

^All  experimental 

data  are  from  Ref. 

7  unless  otherwise 

noted. 

^Averaged  over  the  two  bonds. 


^Assumed  value  in  Ref.  7.  The  uncertainty  is  assigned  in  the 
present  work. 

(iThe  uncertainty  is  assigned  in  the  present  work. 


fiRef.  8. 


V.  CONCLUSIUi^S 


An  HF-SCF-lCAO-MG  study  of  several  acyclical  nitramine  compounds  has 
been  performed  in  the  MNOO  and  MINOO/3  approximations.  The  geometries 
predicted  by  MNOO  tend  to  be  better  than  those  predicted  by  MINDO/3  for  the 
pyramidal  molecules,  but  this  is  not  the  case  for  planar  molecules.  The 
heats  of  formation  predicted  by  MINDO/3  are  in  better  agreement  with 
experiment  than  those  of  MNOO.  Application  of  both  models  to  chlorinated, 
acyclical  nitramines  may  be  undertaken  with  caution. 

Because  of  the  success  of  MINDO/3  in  accounting  for  the  properties  of 
the  model  compound  OMNA,  study  of  the  reaction 
(CH2)N(N02)  -»■  N2OfH2C0, 

which  occurs  in  the  decomposition  of  that  model,  may  be  undertaken  with 
confidence  in  the  results.  Such  calculations  may  be  expected  to  lead  to 
increased  understanding  of  the  dynamics  of  energetic  materials. 
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MINDU/3  fail  to  rotate  the  CNC  plane  with  respect  to  the  N(NU2)  plane. 

This  result  is  consistent  with  the  tendency,  observed  in  this  work,  of 
KINDU/3  to  predict  planar  nitramines.  The  CC^  bond  length  predicted  by 
NNDO  is  in  fair  agreement  with  experiment  while  that  predicted  by  MINOO/3  is 
rather  large.  The  heats  of  formation  predicted  by  both  hNDU  and  mINDO/3  for 
CjiMMNA  are  both  positive  with  the  value  yielded  by  the  former  similar  to  the 
otners  obtained  in  this  work.  It  is  the  first  such  value  obtained  from 
MINDO/3  and  suggests  caution  in  applying  the  model  to  chlorinated  nitramines 

F.  N-chloro-N-methylnitramine  -  Cji(CH2)N(N02) 

Figure  8  contains  a  wedge  diagram  for  N-chloro-N-methylnitramine 
(CiiMNA),  and  Table  6  contains  the  comparison  of  experimental  data  with  MNDO 
results.  The  MINDO/3  calculations  for  CiiMNA  led  to  the  dissociation  of  the 
system  into  NO2  and  (CH2)NCt  fragments.  The  difficulty  may  lie  with 
the  Cl  parameterization  within  MINOO/3.  The  structure  predicted  by  MNDO,  on 
the  other  hand,  is  in  quite  reasonable  agreement  with  experiment.  It  is 
encouraging  to  observe  that  the  CtN  bond  and  CtNN  angle  agree  well  with 
experiment.  The  torsions,  z,CNN0  and  Z.CS.NN0  are  respectively  smaller  and 
larger  than  those  observed.  The  calculated  heat  of  formation  is  positive  in 
agreement  with  the  other  MNDO  values  obtained  in  this  study. 


suggests  the  speculation  that  the  value  is  negative  which  would  be  in 
agreement  with  the  MINDO/3  prediction. 

D.  N-nitro-N-methyl  nitramine  -  IN02){CH2)N(N02) 

Figure  6  contains  a  wedge  oiagram  illustrating  the  structure  of 

N-nitro-N-methyl  nitramine  (oinitromethylamine;  NMNA),  and  Table  4  contains 

a  comparison  of  the  results  of  the  calculations  with  experiment. 

1  2 

Experimentally,  the  two  nitro  groups  and  the  N  N  bond  lie  within  a 

single  plane,  and  the  CN  bond  is  displaced  out  of  the  plane.  Both  MNDO  and 

MINDu/3  yield  the  result  that  the  oxygen  atoms  are  rotated  out  of  the 
12  3 

N  N  N  plane.  The  pyramidal  location  of  the  CN  bond  with  respect  to 
the  plane  of  the  nitrogen  atoms  is  better  reproduced  by  MNDO  than  by 
MINDO/3.  In  the  present  case,  the  NN  bond  length  is  more  than  0.2i  longer 
than  in  nitramide,  a  trend  which  is  more  closely  followed  by  MNDO  than  by 
MINOO/3.  The  CN  bond  lengths  appear  to  be  insensitive,  both  experimentally 
and  theoretically,  to  the  other  moiety  bonded  to  the  amino  nitrogen  atom. 

As  in  the  previous  cases,  MNDO  yields  a  positive  heat  of  formation  while 
MINDO/3  yields  one  that  is  negative. 

E.  Chloromethyl -methyl nitramine  -  (CHgCt) (CH2)N(N02) 

Figure  7  contains  a  wedge  diagram  for  chloromethyl-methylnitramine 
(CjIMMNA),  and  Table  5  contains  the  comparison  between  experiment  and  the 
results  of  the  calculations.  The  bonds  emanating  from  the  amino  nitrogen 
atom  are  rotated  from  the  NINO^)  plane.  In  Table  5,  the  angle  0  is 

c 

defined  to  be  the  angle  between  the  CNC  and  C  NN  planes.  The  MNDO 
calculation  yields  a  result  for  this  quantity  that  is  in  poor  agreement  with 
experiment,  while  MINDO/3  yields  satisfactory  cigreement.  both  MNDO  and 


The 


Q  IQ 

the  structural  data  have  been  obtained  from  electron  diffraction.  * 

third  and  fourth  angles  listed  in  the  table  are  torsion  angles  defined 

according  the  lUPAC  convention  ,  while  the  fifth  angle  may  be  considered 

2  3  4  1 

to  be  the  torsion  for  the  NO  bond  ana  the  0  N  “pseudo-bond".  The 
MNDO  results  are  slightly  better  than  those  obtained  for  nitramide,  while 
those  obtained  from  MINDO/3  are  quite  good.  As  before,  the  NN  bond  length 
yielded  by  MINDO/3  is  appreciably  shorter  than  that  observed  and  may  be 
understood  in  terms  of  the  previously  discussed  population  analysis.  The 
dipole  moment  predicted  by  MNDO  is  in  better  agreement  with  experiment  than 
that  predicted  by  MINDO/3  principally  as  a  result  of  the  latter  predicting 
larger  hybridization  moment  than  the  former.  The  heat  of  formation 
predicted  by  MNDO  is  40  kcal  m"^  too  large  and  of  the  wrong  sign,  while 
that  predicted  by  MINDO/3  is  but  7  kcal  m"^  too  large  and  of  the  correct 
sign. 


C.  Methyl nitr amine  -  H(CH2)N(N02) 

There  is  some  uncertainty  about  the  location  of  the  isolated  H  atom  in 
14 

methylnitramine  (MNA)  ,  and  Fig.  5  contains  a  Newman  projection  of  the 
structure.  The  data  and  the  results  of  the  calculations  are  contained  in 
Table  3.  The  NN  bond  length  is  approximately  the  same  as  in  DMNA  but  0.05A 
shorter  than  in  nitramide.  This  trend  is  reproduced  reasonably  well  by 
MNDO,  but  the  NN  bond  length  predicted  by  MINDO/3  is  too  short  and  does  not 
exhibit  this  behavior.  The  CNNO  torsion  predicted  by  MNDO  is  larger  than 
that  observed  while  MINDO/3  yields  a  planar  structure  for  the  heavy  atoms. 
Although  the  experimental  value  of  the  heat  of  formation  of  MNA  is  not 


known,  consideration  of  (CH2)N02,  H2N(N02),  and  (CH2)2M{N02) 


The  two  angles  provide  a  measure  of  the  agreement  of  the  calculated  and 
observed  descriptions  of  the  pyramidal  structure  of  the  molecule.  The 
alternative  of  constraining  the  heavy  atoms  to  be  coplanar  was  rejected 
because  this  property  should  be  a  consequence  of  the  theory. 

The  results  of  Table  1  indicate  that  MNDO  provides  a  quite  reasonable 
description  of  the  geometry  of  nitramide,  while  that  predicted  by  MINDO/3  is 
rather  poor.  The  NN  bond  length  predicted  by  MINDO/3  is  much  too  short  and 
the  molecule  is  predicted  to  be  planar.  The  short  NN  bond  length  can  be 
understood  in  terms  of  the  HF  wave  function.  The  MINDO/3  Hamiltonian 
results  in  there  being  a  larger  positive  charge  on  the  nitro-N  atom  and  a 
smaller  (larger  magnitude)  negative  charge  on  the  amino-N  atom.  The 
increased  electrostatic  attraction  between  these  atoms  leads  to  a  shortening 
of  the  bond  length. 

The  heat  of  formation  predicted  by  MNDO  is  33  kcal  m*^  too  large  while 
tnat  predicted  by  MINDO/3  is  16  kcal  m“^  too  small  but  of  the  correct 
sign.  It  is  interesting  to  note  that  the  total  molecular  energy  of  MNDO  is 
deeper  than  that  of  MINDO/3.  The  differences  between  the  heats  of  formation 
predicted  by  the  two  models  is  principally  due  to  the  electronic  energies  of 
the  constituent  atoms. 

B.  Dimethylnitramine  -  (CH2)2N(N02) 

Table  2  contains  a  comparison  of  the  results  of  the  calculations  and 
experiment  for  dimethylnitramine  (DMNA).  Figure  4  contains  a  schematic 
illustration  of  the  heavy  atoms *in  the  molecule  which  is  planar.  Most  of 


where  are  the  theoretical  and  experimental  quantities  of  interest 
ana  is  the  uncertainty  of  the  observed  ;  the  total  number  of 
compared  quantities  for  a  molecule  is  N.  The  quantity  differs  from  the 
frequently  encountered  statistical  quantity  y}.  i^o  parameters  are  being 
aajusted  in  oraer  to  fit  experiment,  ana  the  sum  of  the  weights  is  not 
normalized  to  unity.  Excellent  agreement  between  theory  and  experiment  for 
molecular  geometries  would  be  to  within  O.OlX  for  bond  lengths  and  1°  for 
bona  angles. 


A.  Nitramide  -  H2N(N02) 


Structural  information  about  the  most  elementary  nitramine  compound, 
H2N(N02),  has  been  obtained  from  x-ray  scattering^  and  microwave 
spectroscopy,^  The  molecule  is  pyramidal  and  schematically  illustrated  in 
Fig.  3.  Table  1  contains  experimental  data  for  the  molecule  and  the  results 
of  calculations  performed  with  the  MINDO/3  and  MNDO  Hamiltonians. 


The  angle  *  is  definea  experimentally  to  be  the  angle  between  the  HNH 
and  N(NU2)  planes.  It  is  not  theoretically  necessary  for  four  points  to 
lie  in  a  plane,  and  this  possibility  was  allowed  in  the  geometry 
optimizations.  Accordingly,  the  definition 

5(N^N^)X^(N^oS]  *  [t(N^H^)XR(N^H^) 


|$(N^N^)xt(N^oS  II  I 

has  been  adopted.  A  second  angle  'i'  has  been  defined  by  the  relation 
[R(N^O^)}®(N^oh]  •  [R(N^N^)XR(Ahj 

C.O S  'V  *  .  —  -  —  ■  ■  I  -  Ml ■  ■  ■  I  ■  ■  I  III 

|R(N^O^)}a(N^oh  I  I  R(N^N^)XR(N^O^) 
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features  of  electron  correlation  into  a  theory  otherwise  lacking  tne  effects 
of  such  correlation. 


There  are  interesting  aspects  of  the  calculation  of  the  theoretical  heat 

of  formation  of  a  compound  definec  by 

=  U  -  ,  (20 

m  a  01  a 

where  the  total  molecular  energy  is  u  according  to  Eq.  (2),  the  atomic 


electronic  energy  for  isolated  atom  a  is  ,  and  ah^  is  the  heat  of 


formation  of  free  atom  a  from  the  element  in  its  standard  state.  In  both 
hlNDO/S  and  MNOG  the  ah“  are  experimental  quantities  but  the  E®^  are 
theoretical  quantities  calculated  with  two  different  sets  of 
approximations.  It  is,  therefore,  possible  to  perform  calculations  in  the 
two  models  which  agree  very  well  with  each  other  for  Ujj  but  which  yield 
quite  different  values  for  the  experimentally  accessible  heat  of  formation. 

It  is  also  possible  that  ah“°^  is  positive  for  a  stable  compound  which 
illustrates  that  it  should  not  be  confused  with  the  binding  energy 
encountered  in  nuclear  chemistry 

(2i: 


gnuc^^nuc 


In  £q.  (21),  the  total  nuclear  energy  is  for  the  nucleus  consisting  of 
z  protons  of  the  mass  and  N  neutrons  of  mass  M^.  The  binding  energy 
is  always  negative  for  a  stable  system. 


IV.  RESULTS 


The  quality  of  the  theoretical  results  when  compared  to  experiment  is 


measured  by 


electronic  charge  oistribution  to  S.  The  two-center  one-electron  core 
resonance  integrals  are  assumeo  to  be  given  by 

S^j.(ep“  +  s/)/2,  pect  ;  reg  ^8 

c(  3 

where  Sp  and  are  adjustable  parameters  characteristic  of  the  orbitals  and 
their  parent  atoms.  The  two-center  two-electron  repulsion  integrals  are 
determined  by  a  semi-empirical  form  based  upon  a  multipole  expansion  of  the 
integral  in  question.  For  first  row  atoms  calculated  with  a  minimal  basis 
set,  this  implies  that  quadrupole-quadrupole  interactions  are  the  highest 
order  encountered.  This  represents  an  improvement  over  MINDO/3  because  the 
integrals  now  depend  not  only  on  the  atom  pair  but  also  on  the  electronic 
angular  momenta.  The  forms  chosen,  of  course,  yield  the  correct  results  for 
very  large  and  for  very  small  separations  .  Finally,  the  core-core 
repulsion  integrals  are  determined  from  the  quantum  mechanical  monopole 
charge  distributions  of  atoms  «  and  8  modified  by  empirically  determined 
exponential  functions. 

The  semi-empirical  methods  are  parameterized  by  calculating  the  heats  of 
formation,  equilibrium  geometries,  and  electrostatic  dipole  moments  of  a 
number  of  sagaciously  chosen  molecules.  These  calculations  are  repeated 
many  times  in  the  process  of  performing  a  non-linear  least  squares 
optimization  of 

j^Y^Ccalc.)  -Yj^(obs.)j^  wj  ,  (19 

2 

where  the  Y^  are  the  quantities  in  question,  and  the  are  weights 
ascribeo  to  the  Y„ .  This  procedure  creates  an  inter-oependency  among  the 
parameters  which,  because  they  are  empirically  determined,  incorporate  some 
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